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Abstract. OAO 1657—415 is an accreting pulsar in an eclipsing binary system. We analyzed the INTEGRAL core 
program observations of this object and obtained the eclipse light curve in the soft gamma-ray band between 15 
and 40 keV. We note that the gamma rays from the pulsar allow to probe the density profile of the outer layers 
of the B supergiant companion. We find that the density profile of the outer layer can be described by a power 
law with the index a = 8.5. We also note that the fit hints toward smaller inclinations of the system within the 
allowed range 60° < i < 90° llChakrabartv et alJll99.8ll . 
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1. Introduction 


2. Data acquisition and processing 


The accreting pulsar OAO 1657— 415 has been discov ¬ 
ered with the Copernicu s satellite bv IPolidan et al.l ( 1978h . 
IWhite fc Pravdol lll979^ found the 38.22 s pulsations from 
the source. The pulsar spin has been measured by a 
number of observatories and the spin history shows 
large spin-up and spin-down episod es on a timescale 
betw e en a few months a nd a year I Chakrabartv et al.l 
Il993t iBildsten' et al ]|I993). The companion has initially 
been identified as V861 Sco, ye t this identification has 
been proved to be incorrect ifw hite fc Pravdol Il979t 
lArmstrong et alJll98(ilBvrne et al.lll98lll . The pulsar has 
been found to be a member of rare cl ass of eclipsing ob¬ 
jects with an orbital period of 10.4 d llChakrabartv et all 
Il99.3h . and the optical companion identification was made 
possible only recently by the ex cellent X-Ray counterpar t 
location provided by Chandra llChakrabartv et aLll2002ll . 
This identification lead to estimate of the distance at 6.4± 
1.5 kpc, much lower than the llkpc lower limit inferred 
previously with the assumption that the neutron star is 
rotating at the equilibrium period. The low energy X-ray 
spectr um of OAO 1657—415 is very absorbed, N h ~ 10^^ 
cm“^ llPolidan et al.l[l975 iKamata et al.lll990^ . because 
the source lies at a low galactic latitude, and Ginga ob¬ 
servations found a soft exces s below 3JceV_||Kamataet_^D 
119901 Iwhite fc Pravdolll97^ . while iKamata et al.1 lll990h 
found a fluorescence line a 6.6 keV. The hard X-ray spec¬ 
trum can be modeled by a power law steepening at en¬ 
ergies ab ove 20keV, with a hi nt of a cyclotron line at 
« 36keV llOrlandini et al.l[i999ll . 


INTEGRAL is a Gamma-ray satellite equipped with three 
high energy scientific instruments: SPI - gamma spec- 
trom eter, IBIS - gamma imager and JEM-X - X-ray mon¬ 
itor I Winkler et alJl200^ . ISGRI is the low energy detec¬ 
tor on board the IBIS/INTEGRAL telescope - it is a po¬ 
sition sensitive 0.5 meter b y 0.5 meter detector made of 
CdTe llUbertini et al.ll2003l . It works in the energy range 
15 — 1000 keV. ISGRI consists of about 16 384 single de¬ 
tector pixels. The coded mask aperture is placed about 
3 meters above the ISGRI detector. The mask is made 
of the Tungsten segments: about 1 by 1 cm and 1.6 cm 
thick, in 95 by 95 segments system (half open), with the 
basic pattern of 53 by 53. This coded aperture mask tele¬ 
scope collects gamma rays with angular resolution better 
than 12 arcmin and with energy resolution about 10% at 
100 keV with about 2500 cm^ collecting area. 


All the analyzed data come from ISGRI/IBIS detec¬ 
tion layer. We have selected the data intervals from the 
Galactic Plane Scan in the INTEGRAL core programme 
for which OAO 1657—415 was in the fully coded field of 
view, i.e. less than 10° from the axis. 


This condition was satisfied by the observations from 
the the following INTEGRAL revolutions: 44-47, 49- 
56, 58-61, 63 (Spring 2003), and 100, 101, 103, 116, 
118, 119 (Autumn 2003) The total time of data acqui¬ 
sition was ~ 900 ks. The data has been processed us¬ 
ing standard INTEGRAL software provided by ISDC 
(http: //isdc . unige . ch) in version 3.0 with normal con- 
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Fig. 1. The map of OAO 1657—415 region. The source is 
detected at the level of 122 a for a total exposure of 24.5 ks 

figuration for the ISGRI pipeline. A preliminary ana lysis 
of the data has been presented in lPenis et al.l (|20^). 

Since we have used the data segments corresponding 
to different locations on the detector we had to model the 
off axis response of the IBIS/ISGRI detector. We have 
used the Crab observations during the INTEGRAL rev¬ 
olution 170 as our reference. We used two models of the 
detector angular response function: a wide area (± 13° 
square around telescope axis) global fit with a paraboloid 
summed with a two-dimensional Gaussian, and local flat 
surface triangulation. Both gave similar results. 

We have analyzed 480 science windows (INTEGRAL 
science window is a single exposure lasting of 20 min¬ 
utes) and the typical detection significance was 25 sigma. 
We have also detected the pulsations with the period of 
37.16 s. We present the skymap of OAO 1657—415 region 
in Figure^ The source is detected at the level of 122 cr in 
the 15—40keV ene rgy band for a total exposure of 24.5 ks 
llPenis et al ] |2004 . 

3. Eclipse modeling 

The brightness of the source varies from 200-300 mCrab 
to nearly 0 between different observations. We have folded 
the observations with th e orbital period of t he pu lsar of 
10.44days discovered bv IChakrabartv et alJ lll993ll . This 
revealed a beautiful eclipse profile, as can be seen in 
Figure|21 Epoch of the eclipse centre is JP 2 452 705.9. 

The eclipse in the lightcurve is clearly visible. What 
is interesting is that the eclipse ingress and regress are 
gradual, and last almost a day, with the totality, ingress 
and regress, of almost two days. A pulsar is essentially a 
point source in comparison with the B supergiant com¬ 
panion with a radius of « 30 Rq . Such a gradual eclipse 
is expected when the sizes of the sources are comparable, 
or if the eclipsing object is diffuse and has a low optical 



Fig. 2. OAO lightcurve with the best fit of the eclipse 
model (see text for details) 

thickness to the gamma rays of the eclipsed source. In the 
case of OAO 1657—415 light curve we have to consider 
this latter case. 

To model the eclipse we have assumed that the density 
profile of the B supergiant companion can be described as 
a power law: 

(site) ■ 

For a given orbital phase (f> and assumed inclination i we 
can calculate the the column density through the eclipsing 
star by integrating the density profile: 

poo 

ne{4>) = / dr'p{y^r±{(j)f + r'^) (2) 

where ry {(j)) is the projection of distance between the pul¬ 
sar and the companion along the line of sight, and r± (cj)) is 
the projection of this distance on the plane perpendicular 
to the line of sight. At the energies 15 — 40keV the main 
opacity through the star is due to electron scattering. We 
treat the radiative transfer approximately and we model 
the eclipsed flux from the pulsar as 

F{4>) = Fo exp(-ne(0 - (j)o)(TT) (3) 

where Fq is the flux from the pulsar, (j)Q is the center of the 
eclipse and ax is the Thompson cross section. Altough the 
source is intrinsically variable (20-40%), we put Fq con¬ 
stant because the modeled eclipse light curve is a super¬ 
position of many data sets covering randomly the whole 
range of the orbital phase. Such an approach should not 
have an important impact on the fitted values of the 
eclipse parameters (i, 4>o, po, and a), but can increase the 
incertitude on Fq. In the model we use the orbital elements 
inferr ed earlier from the BATSE data llChakrabartv et ^ 
[l99l. The model has five free parameters: i, 4>o, Fq, pq, 
and a. 

The data analysis yielded the statistical errors. We es¬ 
timate the additional systematic errors by requiring that 
the data in the total eclipse part (phases from 0.13 to 
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Fig. 3. The gray contours enclose the regions containing 
68, 95, and 99.9 percent probability. The index a increases 
with decreasing inclination angle 


Fig. 4. The gray contours enclose the regions containing 
68, 95, and 99.9 percent probability. The density po in¬ 
creases with decreasing inclination angle 


0.30) are consistent with zero flux. To this goal we add in 
quadrature a systematic error Osys = 15.5mCrab. This is 
a minimal value which leads /DOF = 1 in this phase 
interval where the flux is zero during the total eclipse. 

To find the best fit model we use the minimaliza- 
tion. We have scanned the parameter space, with the re- 
strictio n 60° < i < 90°, as i nferre d from the timing anal¬ 
ysis bv IChakrabartv et alJ lll993ll . A best fit model light 
curve is presented as the solid line in Figure The model 
curve correspond to the model with the following param¬ 
eters: i = 60°, Fq = 147mCrab, a = 8.6, (j)o = 0.218, and 
logpo = 12.34, where po is in units of cm“^. 

In order to find the confidence regions on the param¬ 
eters we assume that the probability density in the pa¬ 
rameter space is proportional to oc exp(—where 
Pi are the parameters. Such probability density is nor¬ 
malized. We marginalize over the remaining parameters 
to find the probability distributions for each parameter, 
or two dimensional maps of probability density. The best 
fit values with the la conhdence contours obtained using 
the above procedure are a = 8.5±);5, logpo = 12.31±o;3o, 
the center of the eclipse is (po = 0.218 ± 0.001. The for¬ 
mal best fit value for the flux Fq is Fq = 147 ± 1 mCrab, 
and for the inclination i = 60° ±q. The best fit values of 
the parameters corresponds to the peak of the probabil¬ 
ity distribution and the Icr describe the region containing 
68% probability around this peak. The formal error on in¬ 
ferred value of the inclination angle reflects the fact that 
we considered only the angles above 60° and the best fit 
corresponded to the lowest angle considered. 

The relatively large uncertainties in determination of 
the mass index a and the fiducial density pq are due to the 
fact that these quantities are correlated. In FiguresIHland 
0 present the regions containing 68, 95, and 99.9 percent 


of the probability in the space spanned by the index a and 
the inclination angle (Figure|2l, and by the density pq and 
the inclination angle (Figure 0J. Both the mass index a 
and the density pq are anticorrelated with the inclination 
angle of the binary. The increase of the density parame¬ 
ter is due to the fact that as we decrease the inclination 
angle the eclipse probes farther radii from the stellar cen¬ 
ter, while Pq is defined as density at a constant radius of 
30i?Q. The index a increases since at higher inclinations 
the neutron star is eclipsed along a chord line and the den¬ 
sity difference leading to the observed eclipse corresponds 
to smaller difference of the radii from the stellar center. 
Therefore also the index a is correlated with the density 
parameter pQ. 

4. Summary and discussion 

The INTEGRAL observations of the eclipse are the first 
result where a stars outer layers can be probed directly 
by the external gamma ray source. The mass enclosed in 
the region that we probe can be estimated by integrat¬ 
ing equation ^ starting from the point where the eclipse 
becomes total. The result depends on the assumed in¬ 
clination and varies from 5 x 10“^ Mq for i = 90° to 
1.6 X 10“^ Mq for i = 60°. The thickness of the layer 
that is probed with this method also depends on the as¬ 
sumed inclination. Assuming that the width of the total 
eclipse is 0.14 we estimate that with the inclination of 90° 
we probe down to the radius of ^ 20 Rq from the cen¬ 
ter, while at the assumed inclination i = 60° we see the 
structure above 28 Rq. 

The outer layer of the star is the location of the accel¬ 
eration of the wind. Our results can be used along with 
the continuity equation to estimate the wind velocity pro- 
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file: v{r) oc This estimate uses an assumption that 

we already see the accelerating layer. In reality the wind 
may be taking off in the upper portion of the atmosphere 
which is already optically thin to gamma rays. 

Finally we note that the flux outside the eclipse seems 
to increase by a factor of w 2 before entering the eclipse, 
i.e. between the phase 0.4-0.6 and 0.8-1.0. While the data 
looks particularly noisy it seems interesting to investigate 
the possibility that the variation of the flux is due to the 
structure of the outgoing wind. The accreting neutron star 
may leave a hollow trace in the wind and such a trace 
could correspond to decrease of density during a particular 
phase, when the pulsar is viewed through hole in the wind 
material. 

Summarizing, we have presented the first well observed 
case of a gamma-ray eclipse. We show that the eclipse 
lightcurve allows to study the structure of the outer lay¬ 
ers of the companion. INTEGRAL is particularly well pre¬ 
pared for such studies as it allows long duration monitor¬ 
ing of sources. 
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